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Abstract 


This paper provides two types of control relevant models of planer solid oxide fuel cell system with different details. Dynamic models of system 
components which include heat exchanger, reformer and after-burner are also provided along with the necessary formulation of a fuel cell connected 
in parallel with a capacitor. Steady-state and dynamic simulations of fuel cell system for both types of models are performed. The results indicate 
that both models are comparable in predicting stack voltage at lower current load. But, the discrepancy in the stack voltage, power and temperature 


of different components become more prominent at higher current load. 


© 2006 Elsevier B.V. All rights reserved. 
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1. Introduction 


Fuel cells are electrochemical devices that generate electrical 
energy directly from chemical reactions. Unlike heat engine or 
gas turbine where there are losses due to Carnot cycle efficiency 
limitation and mechanical losses, chemical energy is directly 
converted to electrical energy and thus the fuel cell plant effi- 
ciency can be as high as 40-55%. In addition, the heat gener- 
ation from the electrochemical reactions in a fuel cell can also 
be used for cogeneration applications increasing the efficiency 
up to 70%. Some other potential features of fuel cell include no 
moving parts, quiet operation, flexibility in fuel types and size 
of the fuel cell, and most of all rapid load following capability 
[23]. Another attractive feature of fuel cell is that it demonstrates 
same efficiency irrespective of its size, which makes it useful for 
a versatile range of applications including power generation for 
home or a commercial building, mobile equipments and auto- 
mobile. 

A large amount of work has been conducted on the modeling 
of solid oxide fuel cell from the viewpoint of fuel cell design, 
operability and performance, material selection and controller 
design [11,15,24,18,22,8]. The models can also be developed to 
simulate steady-state or dynamic behavior and can range from 
zero to three-dimensional [3,4]. For example, Padulles et al. 
[19] developed a simple model of a fuel cell-based power plant 
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which included species dynamics but no temperature dynamics 
whereas, Hall and Colclaser [10] developed a more detailed 
model for simulating the transient operation of a tubular solid 
oxide fuel cell that included electrochemical, thermal and mass 
flow elements. The latter paper mainly focused on the effect 
of load change on the terminal voltage. The work presented in 
[1,5,15,4,2,12] investigated transient behavior of a planer SOFC 
due to load change using a multi-dimensional, time-dependent 
model and is useful for simulating fuel cell behavior. The mod- 
els are expressed by a set of PDEs which are solved numerically 
to evaluate the performance of the fuel cell based on different 
design parameters. These models, since expressed as PDEs, 
are however not convenient for designing controllers. Sedghisi- 
garchi and Feliachi [21] went a step further by considering 
chemical and thermal aspects of chemical reactions and ohmic, 
activation and concentration losses in the fuel cell stack. The au- 
thors however only considered lumped thermal models which is 
a drawback particularly for a hybrid fuel cell system in conjuga- 
tion with combustor, turbine and heat exchangers for generating 
electricity where different flow streams have different temper- 
ature distributions. Several authors also developed models of 
fuel cell system for analyzing fuel cell performance [7,6,17,13]. 
Table 2.1 of [5] summarizes the chronological progression of 
planar SOFC modeling in terms of the underlying objectives. 
In this paper, two sets of dynamic models of different details 
are presented by considering electro-chemical and thermal as- 
pects of fuel cell. Both of the models are represented by a set 
of first-order ordinary differential equations and thus are useful 
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for simulating transient behavior of the fuel cell as well as de- 
signing model-based controllers. In addition, simplified thermal 
dynamic models of fuel cell system components are developed 
and connected to fuel cell to simulate the dynamic behavior of 
fuel cell system. The fuel cell system dynamic model consists of 
only nonlinear ordinary differential equation and thus can be uti- 
lized to design system-wide controller such as nonlinear MPC. 
The model can also be used to estimate unmeasured states using 
state estimators. In addition, the advantage of using capacitor 
connected in parallel with the fuel cell along with its formula- 
tion is also described. 

The paper is organized as follows: first, general principle of 
the fuel cell is described followed by the modeling objective in 
Section 3; next, two sets of models, lumped and detail model, 
are presented in Section 4 and 5 followed by fuel cell simulation 
in Section 8.1. In this section, both dynamic and steady-state 
results are simulated and analyzed. Next the advantage of using 
a capacitor as an auxiliary power source is described along with 
necessary formulations in Section 6. In Sections 7 and 8.3, the 
fuel cell system model is developed and simulated to analyze 
performance of transient system. 


2. General principle 


A typical fuel cell consists of an electrolyte in contact with 
anode and cathode on either side. Hydrogen rich fuel and air 
are continuously fed into the fuel cell for generating electricity. 
The electrolyte acts as a barrier between anode and cathode al- 
lowing only certain types of ions to pass through it. Fuel cells 
are mainly classified depending on the nature of the electrolytes. 
Some of the most common fuel cells are molten carbonate fuel 
cell (MCEFC), proton exchange membrane fuel cell (PEMFC), 
polymer electrolyte fuel cell (PEFC) and solid oxide fuel cell 
(SOFC) [23]. Solid oxide fuel cell, as its name suggests, uses a 
solid metal oxide as electrolyte. The cell operates at high temper- 
ature which makes it particularly useful for internal reforming 
of natural gases. The high temperature exit gases can be used 
further to generate electricity by a turbine. SOFC can be de- 
signed in different shapes and sizes due to the solid nature of the 
electrolyte, for example, planer and tubular. The high tempera- 
ture of the SOFC however imposes a stringent requirements on 
the material selection due to the thermal expansion mismatch 
between anode, electrolyte, cathode and connector materials. 
Thus, it is important to operate SOFC in such a way that the 
stack temperature remains within the design range. 

In a planer SOFC, several cells are stacked and connected in 
series to complete the circuit. At each cell, hydrogen releases 
electrons at the anode surface which travel through the outer 
circuit and combine with oxygen to produce oxide ions. Elec- 
trolyte which acts as a separator between hydrogen and oxygen, 
and thus prevents direct combustion, allows only certain types 
of positive ions to pass through it. SOFC usually uses Y203- 
stabilized ZO2 (YSZ) as electrolyte which allows oxide ion to 
pass through it to reach the anode surface, where the oxide ion 
combines with Ht to form water. The reactions can be summa- 


rized as follows: 


Anode: H —> 2Ht + 2e (1) 
Cathode: 402 + 2e > 0°% (2) 
Overall: Hz + 502 > H20 + heat (3) 


The electrode potential difference, A E that drives the reaction to 
take place and flow of ions through the circuit can be expressed 
by Nernst equation: 


0.5 

RT H2 Po; 

AE = AEp © 
2F 


(4) 


where A Ep [V] is the standard cell potential; py,, po, and pp o 
the partial pressures of hydrogen, oxygen and steam [atm]; T [K] 
the cell electrode temperature at which the reaction is taking 
place; and R [=1 atm kmol”! K-!] and F [=96485.34 s A/mol] 
are the universal gas constant and the Faraday’s constant, respec- 
tively. The Nernst model in Eq. (4) can thus be used to predict 
cell voltage provided standard cell potential, partial pressures 
and electrode temperature are known which can be derived in 
terms of species and energy balances. 


3. Modeling objective 


Knowledge of the transient and steady-state response of the 
solid oxide fuel cell is important for studying fuel cell perfor- 
mance as well as designing the controller. The level of modeling 
sophistication thus depends on the objective and can vary from 
simple linear state-space models to complex 3-D models. The 
models can also specifically be derived to simulate steady-state 
or dynamic behavior. 

Since in this work, the modeling objective is driven by the 
requirement of a control relevant dynamic model that can predict 
the cell terminal voltage and system component temperatures 
as well as design controller with relative ease, certain model 
characteristics have to be identified: 


e The model should comprise a set of linear or, nonlinear ODEs 
to predict transient behavior of the fuel cell. In other words 
the model should be zero-dimensional to avoid distributed 
model comprising of a set of PDEs which makes it difficult 
to design controllers. 

e It should be able to predict all the important variables includ- 
ing terminal voltage, pressure and temperature. 

e Model should be valid for all operating regions, which is the 
main driving factor for building a first principle model rather 
than a data-based model. 


Based on these modeling objectives, two types of models are 
presented in the following sections. First, a completely “lumped 
model” is derived which assumes uniform temperature through- 
out the cell including both the solid phase and the gas phase. 
Even though this assumption is fairly valid at lower current load, 
where the temperatures of electrode, interconnector and unre- 
acted gases do not differ much, at higher current load, this may 
not be the case. Thus, a second model, named, “detail model” 
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is presented which assumes that different components namely 
electrode, interconnector, fuel and air have different temperature 
distributions. Note that the term “detail model” used in this paper 
is relative to the completely lumped model. In addition, lumped 
thermal models for the balance of fuel cell system, e.g., heat ex- 
changers, reformer and burner are presented and then integrated 
with lumped and detailed model, respectively, to simulate dy- 
namic behavior of fuel cell system stack voltage due to load 
changes and other disturbances. 


4. Lumped model 


The lumped model of a stand-alone planer solid oxide fuel 
cell fed with hydrogen and air is developed based on the follow- 
ing assumptions: 


e The gases are ideal. 

e Pressure is constant inside the gas channel. 

e The ratio of interior and exterior pressures is large enough to 
consider that the orifice is choked. 

e Uniform temperature distribution for the entire fuel cell stack. 

e Ideal mixing of gas inside the fuel cell stack. Hence, the exit 
temperature fuel and air are same as the inside temperature. 

e Heat capacities of the fuel and air are negligible. 

e Negligible heat loss to the surroundings. 

e The Nernst equation can be applied for calculating voltage. 


Based on these assumptions an overall material and energy bal- 
ance is performed around the fuel cell stack as shown in Fig. 1(a). 
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Fig. 1. Schematic diagram of: (a) the lumped model and (b) the detail model 


4.1. Species balance 


The ith component material balance for the fuel cell stack 
can be written as [9,21]: 


PV dx; 4 7 r 
RG a TP TA +A 6) 


where P [atm] and T; [K] are the stack pressure and temperature; 
V [m°] the compartment volume; x; the exit molarity of the ith 
component; and a, a and ñ} are the inlet, outlet and reactive 
molar flow rates of the ith component in mol s~!. 

Component balance for hydrogen in terms of its partial pres- 


sure py, can thus be expressed as: 


dpp, 
dt 


R T; „in - out Ad 
= Ven (m, — An, — ^m) (6) 
where Van [m3] is the anode compartment volume; and nit,» net 
and ñ}, are the inlet, outlet and reactive molar flow rates of H2 
[mol s~!]. The reactive and exit flow rates of hydrogen can be 
expressed by [19]: 


fy, = 2K; (7) 


nen = Ky, Puy (8) 
where K, = No/4F; I [A] the stack current; No the number of 
cells associated in series in the stack; and Ky, is the valve molar 
constant for hydrogen. 

Defining tH, = Van/ Kn, RTs, Eq. (6) can then be written as: 


d — 
q’ = 


(a Kp, PH — 2K,I) (9) 
TH, KH, 2 
Since Ky, is a constant, ty, is a function of cell temper- 
ature only and can be expressed as ty, = TH T*/T; where 
TH = TH, |7,=7*. Eq. (9) thus can be rewritten as: 


dpp, Ts 
dt Ky, TH T* 


(iit, — Km Pm — 2K;2) (10) 


Similarly component balances for O2 and H20 lead to the fol- 
lowing set of equations: 


dpo T, i 
T F = 7s (10, ~ Ko, Po, — Kr!) (11) 
2 “O2 
and, 
dpu,0 T; i 
= F aa 7x 4,0 7 Kmopmo + 2Kr1) (12) 
29 “H20 


4.2. Energy balance 


In the lumped stack modeling, it is assumed that there is no 
temperature variation inside the fuel cell which means that all the 
components of the fuel cell—electrode, interconnector and gases 
inside channels posses the same temperature at any instance. It 
is further assumed that the heat capacity of the gases inside the 
channels are negligible compared to the solid components of the 
fuel cell. Then the dynamic model of the cell temperature, Ty, 
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can be found by performing energy balance around the entire 
fuel cell stack: 


z dT. a Tin 
misC ps - = Sar [ Cp,i(T) dT — So age 
ref 


Ts 

x i Cp,i(T) dT — ny, AH? — VsI (13) 
ref 

where m, and C ps are the mass and average specific heat of fuel 

cell materials excluding gases; Cy; the specific heat of ith fuel 

or air gas entering the system; A AY,° the specific heat of reaction 

of Eq. (3); and V; is the stack voltage. 


4.3. Stack voltage 


Applying Nernst’s equation and considering ohmic, activa- 
tion and concentration polarization, the stack voltage can be 
expressed as, 


Vs = Vo — Nohm — Nact — Neon (14) 
where the open circuit voltage, Vo, is 


RTs Pm PO 


Vo = No AE = No | AE In 
0 0 0 0+ FR 


(15) 
PH20 

Here, standard cell potential AEo exhibits linear relationship 
with cell temperature and can be approximated from the exper- 
imental data provided in Table 2-3 of [23]: 


A Eo (V) = 1.2586 — 0.0002527;, (K) (16) 


Ohmic polarization occurs because of resistance to the flow of 
ions through different components of cell materials. The loss 
can be expressed as: 


Nohm = F(T; MI (17) 


The cell resistance, r(7;), however, is a function of cell electrode 
temperature and can be expressed by the second-order Steinhart 
Hart equation: 


r(T;) = ro exp fe (z — al (18) 


where rọ [Q] is the internal resistance at temperature To [K] and 
æ is a constant. 

Activation and concentration polarization can be calculated 
by the following equations [5,21,14]: 


Nat = a + b log I (19) 
and, 

RT; I 
Neon = OF In (: = x) (20) 


where a and bare the Tafel constant and Tafel slope, respectively; 
and J; is the limiting current. 


5. Detail model 


As stated earlier, at higher current load the assumption of uni- 
form temperature distribution throughout the fuel cell in lumped 


model, as described in Section 4, may not be valid. Thus, in this 
section a relatively more detail model is developed where it is 
assumed that different cell components have different tempera- 
ture distributions and thus the model is “detailed” compared to 
the “lumped” model. The major assumptions that differentiate 
it from the lumped model are described below: 


e Different temperatures among electrode, interconnector, fuel 
and air side gases. 

e Ideal mixing of gases inside the fuel cell channels so that fuel 
and air exit temperatures are same as the temperature inside 
the fuel and air channels, respectively. 

e No temperature variation in the axial direction. 


Based on these assumptions material and energy balances are 
performed around electrode, interconnector, fuel gas channel 
and air gas channel of the fuel cell as depicted in Fig. 1(b). 
Since in this model only temperature distributions are assumed 
to be different, species balance around the cell remains same as 
described in Section 4.1, and so is the cell voltage calculation 
in Section 4.3. The stack temperature T, in Eqs. (15), (16), (18) 
and (20) however has to be replaced by electrode temperature, 
Te. 


5.1. Energy balance around electrode 


The electrode control volume consists of the anode, elec- 
trolyte and the cathode. Even though the electro-chemical re- 
actions take place at the anode and cathode near the surfaces 
of the electrolyte, and thus have temperature variations in the 
direction vertical to the surface area, it can be assumed constant 
due to very small thickness of the electrode. The temperature 
variation along the flow direction is also assumed constant, and 
thus the electrode temperature Te dynamics can be expressed as 
a function of diffusive, convective, radiative and reactive heat 
transfer terms (Fig. 1(b)): 


. dr . ; 
peAc AweCpe ae = (Oar = Os) + (Qia = oO 
—(Onre + QiF + Qnae + Ora) — OT — W (21) 


Here, Ac, Awe, pe and C pe are the electro-chemical surface area, 
thickness, density and specific heat of the electrode material, 
respectively. The fuel side diffusion heat terms in Eq. (21) can 
be expressed as: 


. Tp 
Our = iy, | Cp,H)(T) dT (22) 


ref 


Tp 
O%" = nts, | Cpayo(T) dT (23) 


ref 


where iy, = Ny,9 = 2Kr1, and Tp is the exit temperature of 
the depleted fuel. The air side heat transfer terms in and out of 
the electrode through diffusion in Eq. (21) can be expressed as: 


: Ta 
Qia = 10, | Cp,0,(T) dT (24) 


ref 
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o =0 (25) 


where 10, = K;,I, and Ty, is the exit temperature of depleted 
air. The convective heat transfers between electrode and fuel/air 
gases in Eq. (21) can be written as: 


QhFe = hpeAc(Te = Tr) (26) 
Qnae = hacAc(Te — Ta) (27) 


Here, the fuel side and air side heat transfer coefficients Ape 
and fae are also functions of Tp and Ta, respectively, and are 
evaluated at each instance empirically [16,20]. The radiative heat 
transfer terms between electrode and interconnector of fuel side 
and air side in Eq. (21) can be expressed as: 


OAT — T;') 
(1/€a) + /€i) — 1 
oA(Te — T}) 

A/e) + (1/ei)— 1 

where ø = 5.6704 x 1078 W m~? K~4 is the Boltzman con- 

stant; T; [K] the interconnector temperature; and €a, €c and éi 
are the emmisivity constants of anode, cathode and intercon- 
nector materials, respectively. The geometric factor of radiation 


in this case is assumed to be 1. Last, the heat generation and 
shaft work done by the fuel cell can be expressed as: 


Ow (28) 


Ora = (29) 


Q" = -ih AH? and W= VI (30) 


—hpAc(Te+Ti) 


5.3. Energy balance around fuel side 


The dynamic model of exit fuel temperature, Tp can be ex- 
pressed as: 


y Pi Vchannel C9" dTF (ain — ge) 
RTF dt `F F 
+ (aot — Git.) + (Qhre — Onri) (34) 


Here, heat in, o”, and heat out, og, of fuel cell by fuel flow 
can be expressed as: 


, , Tin 
git = dab | " Cp i(T) AT (35) 
Tref 
Tp 
ort = So apt f Cpi(T) dT (36) 
ref 


where Vee is the temperature of fuel entering the SOFC. The 
diffusive and convective heat transfer terms Qip, ae QhFe 
and Qpri in Eq. (34) are expressed in Eqs. (22), (23), (26) and 
(32), respectively. If the energy stored by fuel gas is assumed to 
be negligible and the specific heats of fuel gases are assumed to 
be constant over the operating temperature region, then Eq. (34) 
reduces to an algebraic equation: 


= r 


where 


Tp (37) 


: LaF = se 
Nr = -QË + —hpeAcTe — hpi AiT + |- NO AMC + Ai, oC tno — ^in Cot | Tref 


out 


— àt mout “ro r 
Dp = AH oĆCp, m0 — np, Cpto 


5.2. Energy balance around interconnector 


The interconnector temperature, 7; dynamics can be ex- 
pressed by the following equation: 


dT; 
piAi Aw, C pi rs = Qnri + Ore + Qhai + Ora (31) 


d 
where pi, Aj, Awi, C pi and T; are the density, area, thickness, 
specific heat and temperature of the interconnect materials, re- 
spectively. Convective heat transfers between interconnect and 
fuel/air in Eq. (31) can be expressed as: 


Onri = hri Ai(Te™ — Ti) (32) 
Onai = hai Ai (TR — Ti) (33) 


where hpi and hai are fuel side and air side convective heat 
transfer coefficients adjacent to the interconnector surface, re- 
spectively, which are estimated empirically at each instance. The 
radiative heat transfer terms in Eq. (31) are expressed in Eqs. (28) 
and (29). For planer SOFC the area A; can be approximated as 
the electro-chemical surface area, Ac. 


. out A 
— > ne CH —hpeAc — hpi Aj 
ee 


—2hpAc 
5.4. Energy balance around air side 


The dynamic model of air exit temperature, Ta can be ex- 
pressed as: 


Pi Vehannel Co" dTa . 
2 RT, ar T (28 OR) 
+ (got — O) + (Qnae — Onai) (38) 


Here, on and oy are heat in and out of fuel cell through air 
flow and can be expressed as: 


: f Tin 
gin = ra f * Cy (TAT (39) 
Tref 
out _ „out TA : 
On = Seas Cp,i(T) dT (40) 
Tref 


: : : in out 
The diffusive and convective heat transfer terms Qa» Qaa. 


Qnae and Qna; in Eq. (38) are expressed in Eqs. (24), (25), 
(27) and (33). Assuming energy stored by the air is negligible 
and constant specific heat of air gases for the operating temper- 
ature range, the dynamic equation of Ta can also be reduced to 
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an algebraic equation: 
Na 


Ta = 
A=, 


(41) 
where 


—haAc(Te+T;) 


E = —Tw——_—, 
Na = -O8 - b3 APON — AG, EM, | Tret -hacAcTe — haiAih 


„out mout „r pout 
DA = j hi Chi = 70) Cp, P Ac Ac = haji Aj 
Sx  —— 


—2hA Ac 
6. SOFC along with capacitor 


For a stand-alone solid oxide fuel cell, the model equations 
describing the fuel cell behavior are expressed by the species 
balance [Eqs. (10)—(12)], energy balance [Eqs. (13), (21), (31), 
(34) and (38)] and the stack voltage equation [Eq. (14)]. These 
equations can be put in a concise form by: 


x= f(x, i, Ta, Vs) (42) 
V; = Vo(x) — r(x)la (43) 


where x contains the states describing the partial pressures and 
temperatures of the SOFC models; ñ; the inlet flow rates of fuel 
or air; I4 the demand current load which passes through the fuel 
cell; and V; is the stack voltage produced by the fuel cell. For 
simplicity, only ohmic loss has been taken into consideration in 
the above output equation. For a stand-alone fuel cell, the fuel 
cell current Ifc equals the demand current Jy. From the output 
equation, it is evident that a sudden change in the demand cur- 
rent will be associated with an instantaneous change in the stack 
voltage as shown by the transient responses (Fig. 6). This in- 
stantaneous change in the stack voltage can not be avoided no 
matter what type of advanced control is used due to the con- 
straints on the manipulated variables that are typically the inlet 
fuel and air flow rates. To avoid this sudden loss in voltage and 
possible damage to electrical equipment, an ultra-capacitor of 
sufficient capacity can be used in parallel with the fuel cell as 
an auxiliary power source as shown in Fig. 2. The advantage of 
the capacitor can be seen intuitively—when there is a sudden 
change in the demand current, the capacitor will share the load 
and provide additional power. Thus, instead of sudden drop in 
the stack voltage, it drops smoothly depending on the capac- 
itance of the ultra-capacitor. This gives an added boost to the 
controller connected to the SOFC system to keep the voltage at 
its referenced value. By avoiding sudden drop of the voltage the 
controller copes with only the slow change of the voltage and can 
bring the voltage at its reference value by increasing fuel flow 
rates within its constraints more easily. Mathematically this can 
be expressed as follows. 

Assume that the current through the capacitor and the fuel 
cell are Jy, and Ife, respectively. Then Ife can be expressed in 
terms of the demand current as, 


Ir = Ia — Tac (44) 


From the definition of capacitance, the current passing through 
an ultra-capacitor connected in parallel with a dc voltage source 


Fig. 2. Fuel cell along with a capacitor. 


can be expressed as, 


dV, 
Iu = —C a (45) 
where C is the capacitance [F] and V; is the terminal voltage 
across the capacitor and the dc voltage source. Then current 
through the fuel cell can be rewritten as, 


dV; 
lig = Ig + C L (46) 


Stack voltage of the fuel cell as described by the output equation 
(43) then can be written as, 
dV. 
rC —* = VW —rly— Vs (47) 
dt 
From Eq. (47), it is seen that the stack voltage has now taken 
the form of a first-order ODE instead of a static nonlinear output 
equation. Thus, instead of sudden voltage drop for a step increase 
in the demand current, the voltage drop slows down depending 
on the term rC, as shown in Section 8.2. 


7. Fuel cell system 


In the previous sections, a stand-alone fuel cell is described 
and modeled which can be used for simulating open loop tran- 
sient behavior or designing model-based control. The SOFC 
under consideration is fed with H3 as fuel and air as a source 
of O2. But in practical situations a fuel cell is often associated 
with other components so that the SOFC plant can be fed with 
natural gas or other hydrogen rich fuel instead of pure hydro- 
gen. One such SOFC system is shown in Fig. 3 which includes 
a reformer for converting methane (CH4) into hydrogen, two 
heat exchangers for preheating fuel and air before feeding them 
into the reformer and fuel cell stack, respectively, one burner 
for burning unreacted fuel and, compressors for blowing fuel, 
air and steam. Often the fuel cell system includes a turbine to 
generate auxiliary power from the exhaust gas, fuel processing 
unit (e.g., desulphurizer) and a power conditioning unit (e.g., a 
voltage regulator, a dc/dc or dc/ac converter). These components 
are an integral part of the fuel cell system and called as balance 
of plant (BOP). Examples of process flow diagrams of SOFC 
systems are given in Fig. 1-14 of [23] and in Figs. 3.1, 3.2, 3.12 
and 3.14 of [5]. 
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Fig. 3. SOFC system with heat exchangers, reformer, burner and splitter. 


In the SOFC system that is simulated in this work (Fig. 3), 
methane is pressurized and fed to the fuel heat exchanger for 
preheating by the exhaust gas from the burner. The preheated 
methane then enters an external reformer along with steam where 
methane is being converted to hydrogen through reforming and 
water—gas shift reaction. The product gas from reformer en- 
ters into the anode compartment of the fuel cell stack. Pres- 
surized air is also preheated in another heat exchanger by the 
hot exit gas from the fuel heat exchanger and sent to the cath- 
ode compartment of the fuel cell stack. Hydrogen from the 
anode compartment and oxygen from the cathode takes part 
into the electrochemical reactions to produce power at the elec- 
trode. The depleted fuel and air from the fuel cell stack is then 
fed into a burner to produce heat from the unreacted methane, 
hydrogen and carbon monoxide. The exhaust from the burner 
is then sent to the fuel and air heat exchanger consecutively 
as described earlier. The exhaust gas from the air heat ex- 
changer is then sent for heat recovery in the form of steam 
and hot water. The following sections provide a brief descrip- 
tion and thermal model of different components of the fuel cell 
system. 

Once the system model has been developed, it can be used 
to design model-based controller. The designing of the con- 
troller is dependent on the control objective and thus the key 
manipulated and controlled variables can also differ. For exam- 
ple, if the control objective is to improve load following per- 
formance of the system, then the stack voltage can work as a 
controlled variable whereas the fuel, steam and air (or, oxygen) 
flow rate to the fuel cell and the burner acts as manipulated 
variables for the controller. The flow rates can be controlled 
by using appropriate compressors, blowers and valves which 
can also be mathematically represented by ODE models. For 


this particular case, the split ratio of the exhaust gas to the heat 
recovery unit and heat exchangers can also act as a manipu- 
lated variable. The load and flow temperatures act as distur- 
bances for the system. The detailed design of the controllers is 
however out of scope of this paper and thus is not discussed 
herein. 


7.1. Fuel and air heat exchangers 


Two heat exchangers have been used for preheating fuel and 
air in the fuel cell system. Both heat exchangers are assumed 
to be counter-current double pipe heat exchangers. A portion of 
hot stream which is the exhaust gas from burner is fed into the 
outer pipe and the cold stream is fed into the inner tube counter- 
currently. Assuming that the heat exchanger operates at constant 
pressure (a mild assumption), the general thermal model of the 
heat exchanger can be expressed as follows: 

oT oT 
pCypA — = —pCpvA — + DQ (48) 
ot dz 
where p is the gas density, Cp the heat capacity of the gas at 
constant pressure, A the cross-sectional area and Q is the heat 
transfer rate per unit area based on the heat transfer area xD. 

Since the heat exchanger model as expressed in Eq. (48) is 
a partial differential equation, it is not a control relevant model. 
The model is thus converted into a set of ordinary differential 
equations by dividing the heat exchanger into n nodes along the 
flow direction as shown in Fig. 4. It is assumed that each node 
acts as a continuous stirred tank heater so that the temperature 
inside each node is same as the exit temperature of the gas of 
that particular node. Then energy balance of ith node for the tube 
(cold gas) and shell (hot gas) can be written as, 
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Fig. 4. Heat exchanger divided into n nodes along the length. 


Table 1 

Model parameters 

Parameter Value Unit 

Number of cells, No 384 

Cell area, Ae 100 cm? 

Kp, 8.43 x 1074 kmol atm™! s7! 
Kmo 2.81 x 1074 kmol atm! s7! 
Ko, 2.52 x 1073 kmolatm7! s7! 
TH 26.1 s 

TH o 78.3 s 

ca 2.91 s 

T* 1000 aC 

Constant resistance, ro 0.126 Q 

Resistance slope, œ —2870 

Electrode thickness, he 0.25 mm 
Interconnector thickness, hi 1.5 mm 

Electrode density, pe 6.6 gcm? 
Interconnector density, p; 6.11 gem? 
Specific heat, Cpe, Cpi, Cps 0.4 Jg! K! 

Heat of reaction, AÑ; —0.2418 x 10° Jkmol7! 
Emissivity, €a, €c, €i 0.9 

Table 2 


Balance of plant parameters 


Air heat exchanger: 

Di tube = 0.20 m, Do tube = 0.205 m, Di shen = 0.40 m, L = 10m 
Fuel heat exchanger: 

Di tube = 0.05 m, Do,tube = 0.055 m, Di shen = 0.10 m, L = 200m 
Reformer: 

Reaction area, Arx = 1000 m2, reformer volume, Vr = 10 m? 
Burner: 

Burner volume, Vg = 1 m? 


Tube: 
dT.(i) TÐ — T(i — 1) 
PeCpcAcx a — = —PeC pede Mex 
+7 DoUo(Th(i) — Te(i)) (49) 
Shell: 
dT (i) Thi + 1) — Th 
PrC ph Abx — 4 — = — Pn ph Uh Ax a 
— x DoUo(Th(i) — T.(i)) (50) 


Here, subscripts ‘c’ and ‘h’ stand for the cold side and hot side 
energy balance; Acx and Apx for the cold and hot side flow cross- 
sectional area; Az the length of each node; Do the outer diameter 
of the tube; and Uo! is the overall heat transfer coefficient based 
on Do. The total number of ordinary differential equations for 
the lumped heat exchanger model in this case is 2n. 
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Fig. 5. (a) Power—current, (b) voltage—current and (c) temperature—current 
steady-state curves for planer SOFC. 


7.2. Reformer 


Preheated methane from the fuel heat exchanger and a sepa- 
rate stream of steam are fed into the reformer where endothermic 
reaction takes place between CH4 and H20. 

Reforming reaction: 


CH4 + H20 > CO +3H2, AÂÎÂg? = 206.1kJmol! (51) 
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Fig. 6. Transient responses of voltage, power and temperature due to a load change of 500-600 A. Here, ñy, = 5 mol s7!; no, = 10 mol st; Ti and re are 700, 


800 and 1000 °C for (a and b), (c and d) and (e and f), respectively. 


If the reforming reaction rate is expressed by 7p, then material 
balance for the reactant and product components are, 


nen, = ney, — PR (52) 
ALO = hho — FR (53) 
co =FR (54) 
in, = 37R (55) 


where 

; Ea aay 

rR = koPcu, exp | -—— |] Arx mols (56) 
RTR 

ko = 4274 mol m~? s7! bar7! (57) 

Ex = 82 kJ mol! (58) 


The mixture of CO and H20 then take part in the reversible 
exothermic shift reaction to produce more hydrogen. 
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Fig. 7. Transient responses of voltage and temperature due to a change in hy- 
drogen flow rate from 5 to 6 mol s~!. Here, no, = 10 mol s71; Toei and Ti, are 
700°C. 


Water-—gas shift reaction: 
CO+H20 = CO: +H, AAs? = —41.15kJmol7! (59) 


If the extent of reaction and equilibrium constant for the above 
reaction are expressed by & and Ks, respectively, then the mate- 
rial balance for the reactant and product components are, 


Ape = ùm +E (60) 
ùO, = £ (61) 
ho = Îco — Ẹ (62) 
ipo = imo — & (63) 
where 

Ks(Te) = Eln, + &) (64) 


(aco — €)(tH,0 — £) 
The equilibrium constant Ks can be expressed as a function of 


temperature [4], 


5693.5 
Ks(Tp) = exp (= + 1.077 In Tr + 5.44 x 10°-*7p 
R 


49, 170 
TR 


— 1.125 x 10-773 — - 13.148) (65) 
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Fig. 8. Transient responses of voltage and temperature due to a oxygen flow rate 
change of 10 to 12 mol s~!. Here, ny, = 5 mol sl; Tael and Ti are 700°C. 


Now that Ks is known, the extent of reaction can be calculated 
by solving Eq. (64). The energy balance around the reformer 
can then be written as, 
PYRO = ACH, f Cp,cu, (T) dT 

ref 


in 


„in Tho - out TR 
+1ito f Cp,H,o(T) dT — X ny | Cp i(T) dT 
Tref Tref 
+ rrAAr° +é AHs° (66) 


where Tp is the reformer operating temperature. 
7.3. Burner 


The unreacted fuel and air from the fuel cell is then combusted 
together for heat recovery. In this stage it is important to supply 
additional oxygen so that all the unreacted fuel from fuel cell 
can be consumed. This is done by directly feeding the burner 
with another air stream. This also gives an additional degree of 
freedom for the controller to maintain system temperature within 
operating range. The underlying assumptions for building the 
burner model is given below: 


e Enough Oz is supplied so that all depleted gas from the fuel 
cell can be consumed. 
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air 


e Ideal gas mixing inside the burner chamber so that the exit 
temperature of the burner is same as the inside temperature. 
e Burner operates at constant pressure. 


Energy balance around the burner, 
. dTp ee, [E 
pVpCp 77 Ere f Cpi(T) dT 
ref 


i Ta Tg 
+ an J Cp,(T) AT -X i I Cuar 
Tref Tref 


„in A -in A on a 
— ay, AH, m’ — ACH, A H; ca," — ACO A H, co? (67) 


where 7p is the burner temperature; subscripts ‘F’, ‘A’ and ‘B’ 
stand for fuel, air and burner, respectively; the specific heat C p 
and the density p are also calculated empirically at each com- 
putational instance. 


8. Simulation result 


The lumped and detailed models deduced in Sections 4 and 5 
can be used for simulating fuel cell behavior under different op- 
erating conditions. By considering reforming reaction kinetics, 
it can also be extended to include internal reforming and thus be 
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Fig. 10. Transient responses of voltage, power and temperature due to tempera- 
ture change of air from 700 to 900 °C. Here, ñy, = 5 mol svi; no, = 10 mol s7!; 
current load is 500 A; and Te] is Kept constant at 700°C. 


used to simulate SOFC behavior with natural gas as fuel. Since 
these models comprise a set of nonlinear ODEs, they can also 
be used to design model-based controller. 

The objective of this paper is to synthesize the SOFC cell 
model with that of BOPs to form a control relevant dynamic 
model in system level. The synthesis is by no means a simple 
assembly of the existing models. For example, in order for the 
system to be controllable, additional component with the asso- 
ciated differential equation is needed to avoid the sudden drop 
of voltage due to step change of load. In addition, to make the 
complex system model control relevant, many component mod- 
els including the cell model have to be re-established. However, 
the model parameters needed to be populated from experimental 
and/or literature data for simulation and controller design. Most 
of these data are extracted from [19,21,5,23,20]. Some model 
parameters are not available and thus have been estimated em- 
pirically within sensible limits. The parameters of the balance 
of plant components can be estimated based on the capacity of 
the respective units using engineering design practice. In this 
work, the key parameters are estimated iteratively so that the 
operating conditions can be met at steady-state. Tables 1 and 2 
contain the parameters of a planar SOFC stack made from 384 
cells connected in series and the parameters of a SOFC system. 
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8.1. Fuel cell 


The lumped and detailed models are used to generate power— 
current, voltage—current and temperature—current steady-state 
curves for a fixed 85% fuel utilization? and excess oxygen flow 
ratio of four times stoichiometric value. The inlet temperatures 
of fuel and air were kept constant at 700°C. The simulations 
were done using MATLAB 7.0 by considering the ohmic loss in 
order to compare the results with species-only model described 
in [19]. The existing solvers in MATLAB however can not solve 
differential algebraic equation (DAE) of DAE index greater than 
1. Thus, the exit temperatures of fuel and air gases in the detail 
model have been calculated using the ODEs [Eqs. (34) and (38)] 
instead of the algebraic equations expressed by Eqs. (37) and 
(41), respectively. 

Fig. 5(a and b) shows the steady-state power and voltage out- 
put from the lumped, detail and species-only dynamic models. 
Here, the operating temperature of the species-only dynamic 
model [19] is assumed to be 900°C. From the figures, it is ev- 
ident that the power and voltage predictions from the proposed 
models are comparable whereas the species-only dynamic model 
predicts higher cell terminal power and voltage. In fact, the stack 
voltage and power predicted by the species-only model depends 
on the cell operating temperature which is dependent on the fuel 
and air inlet temperature as well as the demand current. Thus, 
instead of being constant, it varies with operating conditions and 
thus the assumption of constant fuel cell temperature is not valid. 

The stack temperature from the lumped model and the tem- 
peratures of electrode, interconnector, depleted fuel and air gases 
are shown in Fig. 5(c). From the figure it is seen that with in- 
crease in current load fuel cell operating temperature increases. 
It is also to be noted that the difference in exit temperatures of 
fuel and air gases increases with increase in current load which 
are different from the uniform cell temperature of the lumped 
model. Since stack voltage is a function of the fuel cell tempera- 
ture, the output stack voltage from the fuel cell described by the 
lumped and detail model also differs at higher demand current 
as shown in Fig. 5(b). 

Simulations were also conducted to compare the transient 
responses of the fuel cell described by the lumped and detail 
models. The transient response of the fuel cell in terms of fuel 
cell stack voltage and cell temperature has been investigated by 
changing the current loads, fuel and air flow rates as well as their 
inlet temperatures (Figs. 6-10). The transient response of fuel 
cell stack voltage and cell temperature expressed by the models 
due to an increase in current load from 500 to 600 A is shown 
in Fig. 6. Here fuel and O2 flow rates are kept constant at 5 and 
10 mols~!, respectively. The inlet temperatures of fuel and air 
however are 700°C for Fig. 6(a and b), 800°C for Fig. 6(c and 
d) and 1000°C for Fig. 6(e and f). 

From the simulation results, it is noticed that the dynamic 
behavior of the fuel cell stack voltage is greatly dependent on 
the inlet temperatures of fuel and air entering the fuel cell. For 
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Fig. 11. (a) Voltage response of SOFC connected in parallel with a capacitor. 
Demand current sharing by (b) the fuel cell and (c) the capacitor. 


all the cases, there is a sudden drop in the voltage associated with 
the step increase in the load. This sudden voltage drop comes 
from the ohmic loss term r(T)/ in the expression of stack voltage 
Eq. (14). However, since the fuel cell stack temperature, Ts, in 
the lumped model and electrode temperature, Tę, in the detailed 
model are function of the current load, they also increase with 
an increase in the load (Fig. 6(b, d and f)). The increase in the 
cell temperature is associated with a decrease of both internal 
resistance, r(7) in the ohmic loss term and the standard electrode 
potential, A Eo as shown in Eqs. (18) and (16), respectively. 
Thus, there is a net gain or loss in the stack voltage depending 
on the magnitude of ohmic loss and standard electrode potential. 
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Fig. 13. Transient responses of fuel cell system due to change in methane flow rate from 3 to 4mols~!. Here, current load 7 = 500 A; no, = 6mol sl; ùm o = 
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For low fuel and air inlet temperatures (700°C), there is a net 
gain of stack voltage from the initial drop (Fig. 6(a)). However, 
as the inlet temperatures of fuel and air are increased, the gain of 
voltage due to decrease in the ohmic resistance is overshadowed 
by the decrease in the standard electrode potential. This behavior 
results in the gradual decreasing of the stack voltage from the 
initial voltage drop as shown in Fig. 6(c and e). 

The transient behavior of the fuel cell in terms of voltage and 
cell temperature for a step change in the fuel flow rate, ñy, , from 
5 to 6 mol s7! is shown in Fig. 7. The current load J, O2 flow rate 
no, and the inlet temperature of fuel and air are kept constant 
at 500 A, 10 mols~! and 700°C, respectively. Figs. 8-10 show 
the transient responses for a step change in O2 flow rate from 10 
to 12mols~!, and step changes in the inlet temperature of fuel 
and air from 700 to 900°C, respectively. 


8.2. Fuel cell with capacitor 


For a stand-alone SOFC, the voltage drops suddenly when 
the load was increased from 500 to 600A [see Fig. 7]. How- 
ever, when the fuel cell is connected in parallel with a capacitor, 
the combined system behaves as a first-order system. As the ca- 
pacitance of the capacitor is increased from 0 to 500F for the 
combined system, the load sharing capability of the capacitor 
increases [Fig. 11(b and c)] and the sharp decline in the voltage 
is smoothed out [Fig. 1 1(a)]. In other words the capacitor works 
like a filter for the demand current. The final steady-state stack 
voltage for all the cases are same which is evident from Eq. (47). 


8.3. Fuel cell system simulation 


The fuel cell system described in the previous section is sim- 
ulated for a demand current change of 500-600 A keeping all 
other input flow rates and disturbances constant. The balance 
of plant parameters, used in the simulation, are designed so 
that operating constraints are met and are presented in Table 
2. The dynamic response of the fuel cell system in terms of 
stack voltage and different component temperatures is shown in 
Fig. 12. Since the stack voltage of the fuel cell system is de- 
pendent on the input flow rates and temperatures, the prediction 
of voltage is not straightforward. For example, for this case in- 
crease in demand current is associated with increase in the fuel 
and air inlet temperature of the fuel cell as well as increase in 
hydrogen production rate and decrease in steam out of the re- 
former. But since partial pressures are not only a function of 
the corresponding flow rates, but function of temperatures as 
well, the partial pressures of H2 and H2O decreased and in- 
creased, respectively, resulting in additional loss in the stack 
voltage. 

For an increase in the methane flow rate, the stack voltage 
from the fuel cell system increases due to increase and decrease 
in the hydrogen and steam partial pressures, respectively, as 
shown in Fig. 13. Here, the hydrogen production from reformer 
increases as well as the reformer temperature, fuel cell inlet flow 
temperatures and burner exhaust temperature. The increase in 
hydrogen flow rate increases the partial pressure which finally 
results in the increase of stack voltage. 


9. Conclusion 


The main objective of this work was to develop control rele- 
vant model of the fuel cell system so that it can be used to design 
model-based controller. For this purpose, two sets of compara- 
ble models, lumped and detail model both represented by zero- 
dimensional nonlinear ODEs, have been developed. Simplified 
thermal models of system components have also been provided. 
The advantage of using capacitor in parallel with the fuel cell is 
discussed with appropriate formulation and simulations. Simu- 
lations are conducted for both of the models and entire system 
to compare steady-state and dynamic behavior. From the simu- 
lations as well as the understanding of the model, the following 
conclusions can be drawn: 


e SOFC temperature is an important variable in predicting stack 
voltage and thus species-only model which assumes constant 
operating temperature is not valid unless an isothermal oper- 
ation is enforced by other heating device. 

e Stack voltage, current load and cell temperature are internally 
dependent. Change of one affects all of the others. This is 
more prominent in the fuel cell system with after-burner and 
heat exchangers. Change in fuel flow rate or demand current 
or even inlet temperature leads to a change in hydrogen flow 
rate in the fuel cell from the reformer. This leads to a change 
in fuel cell temperature and stack voltage. The unreacted fuel 
then goes to after-burner and heat exchangers which again 
affects the temperatures of the reformer and fuel cell. The 
effect of flow rates or disturbances thus propagates through 
the entire fuel cell system. 

e Steady-state stack voltage, power and temperature of the fuel 
cell predicted by both of the models are comparable at lower 
current load. But as the current load is increased they differ 
in magnitude and thus it is important to choose appropriate 
model depending on the control objective. For example, for 
designing a controller where computational power is limited, 
the lumped model can be used. But for evaluating perfor- 
mance of the fuel cell, detail model is preferable. In addition, 
the lumped model considers that the exit temperature of de- 
pleted fuel and air are same. In fact, these temperatures differ 
widely for higher current load and thus detail model should 
be used, especially when the fuel cell is connected to an after- 
burner for recovering energy. 

e At lower operating temperature, the stack voltage initially 
drops suddenly for a step up of the current load and then gains 
some voltage due to decrease in ohmic loss. But at higher op- 
erating temperature, the gain in stack voltage due to lower 
ohmic loss is overshadowed by the decrease in standard elec- 
trode potential resulting in more decrease in the stack voltage 
from the initial drop. Thus, fuel cell may show different dy- 
namic characteristics depending on the operating region. 

e For the detail model exit temperatures of fuel and air can 
be expressed either by ODEs or by DAEs. However, solving 
DAE may require additional computational power and may 
not be helpful if the model is used to compute future control 
actions using optimization algorithms such as MPC. Thus, 
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in our simulation ODEs have been used to simulation exit 
temperature of fuel and air out of fuel cell. 

e Capacitor can act as an auxiliary power source for the fuel cell 
system soothing out sudden change in the voltage. It can help 
in two ways: first, electrical equipments will be less prone 
to damage due to the sudden fluctuation in the voltage, and 
second, any controller responsible for maintaining the voltage 
will be able to bring the voltage to its reference value more 
easily. 

e It is to be noticed that the lumped and detail models are de- 
scribed by nonlinear ODEs with 4 and 7 states, respectively. 
Hence, during simulation the performance of lumped and de- 
tail models do not differ much in terms of computational time. 
For example, the cpu time to simulate the transient response 
for load change described by lumped and detail models are 
0.6409 and 0.8012 s, respectively [Fig. 6(a and b)]. The com- 
putation, in this particular case, is performed using Dell In- 
spiron 600 M with 2.00 GHz Intel Pentium Mobile processor 
and 1 GB RAM running Microsoft Windows XP Home as 
an operating system. For control purpose, which is not pre- 
sented in this paper, the cpu time for linear controller (e.g., 
linear MPC) does not differ much. But for nonlinear con- 
troller obviously the computational time is higher based on 
the sample time to estimate unmeasured states using non- 
linear state estimator (e.g., unscented Kalman filter) and the 
sample time of the controller (e.g., nonlinear MPC). However, 
with the advent of cheap computational power and dedicated 
servers for applying advanced controller, this is not a problem 
anymore. 


Overall, the developed models can predict all the important 
variables of the fuel cell system. In addition the models being 
zero-dimensional in nature can be used for designing model- 
based controller such as linear and nonlinear MPC. Moreover 
since usually all the states are not measurable, the same model 
can be used for estimating states using state estimator such 
as extended Kalman filter (EKF) or, unscented Kalman filter 
(UKF). 
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